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Abstract 

Gamma-ray bursts (GRBs) are the most luminous electromagnetic explosions in the Universe, which emit up to 8.8 x 10 54 


erg 


isotropic equivalent energy in the hard X-ray band. The high luminosity makes them detectable out to the largest distances yet 
explored in the Universe. GRBs, as bright beacons in the deep Universe, would be the ideal tool to probe the properties of 
high-redshift universe: including the cosmic expansion and dark energy, star formation rate, the reionization epoch and the metal 
enrichment history of the Universe. In this article, we review the luminosity correlations of GRBs, and implications for constraining 
the cosmological parameters and dark energy. Observations show that the progenitors of long GRBs are massive stars. So it is 

m expected that long GRBs are tracers of star formation rate. We also review the high-redshift star formation rate derived from GRBs, 
and implications for the cosmic reionization history. The afterglows of GRBs generally have broken power-law spectra, so it is 

w possible to extract intergalactic medium (IGM) absorption features. We also present the capability of high-redshift GRBs to probe 
the pre-galactic metal enrichment and the first stars. 
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1. Introduction 

Gamma-ray bursts (GRBs) are among the most intrig uing 
phenomena in the Universe (for reviews, see Meszaros 2006; 
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IZhand2007 : Gehrels et alJl2009 : Kumar and Zhand2014 ). Ac 

cording to the duration time Tg o, GRBs are usually classified 
into two classes: long GRBs (T^ > 2 s) and short GRBs 


(Tgo < 2 s) (iKouveliotou et alJl!993h . Long GRBs are thought 
to arise when a massive star (> 25M Q ) undergoes core col¬ 
lapse, but the progenitors of short GRBs are mergers of dou- 
ble neutron star or a neu tron star and a black hole binary 
( Wooslev and Bloom 2006 1. Long GRBs can be made “rela¬ 
tive standard candles”, using luminosity correlations that have 


X 


been found in prompt and afterglow phases (i.e.. lAmati et al 
12002 : Ghirlanda et al. 2004 at Liang and Zhandl2005l) . In view 


that the history of the Universe during the so-called “dark 
age” (from cosmic background radiation at z ~ 1100 to the 
epoch when fi rst stars were formed a round z ~ 20) is still 
poorly known ( Barkana and Loebll2001 ). GRBs, as bright bea¬ 
cons in the deep Universe, would be the unique tool to illu¬ 
minate the dark Universe and allow us to unveil the reion¬ 
ization history. GRBs provide ideal probes of the formation 
rate and environmental impact of stars in the high-redshift 
universe, including the metal enrichment of the intergalactic 
medium (IGM). Meanwhile, the infrared (IR) and near-IR af¬ 
terglows of long G RBs are also expected to be detectable out to 
very high redshifts ( LaqibandReichart 20 0(3; Ciardi and Loehl 


brighter. The afterglow intensities also decrease with time. So 
according to the standard theory of GRBs, there is little or no 
decrease in the flux of GRB afterglows at a given observed time 
with increasing redshift for a single GRB. Consequently, GRBs 
can be used as powerful probes of the very high redshift Uni¬ 
verse. The prospects of using GRBs as a cosmological tool are 
exciting: 

• The high luminosities of GRBs make them detectable out 
to high redshifts. 

• Gamma-ray photons suffer from no extinction when they 
propagate towards us. But the optical photons from super¬ 
novae will suffer extinction from the interstellar medium 
(ISM); 

• The correlations between GRB spectral properties and 
energetics have been shown to be powerful tools that 
“standardize” GRB energetics. So GRBs can be used 
to constrain cos mological parameters and the nature of 
dark energy (i.e., Dai_etaljj2004^ Ghirlanda et al. 2004h ; 


Liang and Zhang82005t Ghirlanda et aljl2006a): 


200(3 Bromm and Loeb Eooi : Iwang et al. 20121) . The 

is that the cosmological time dilation translates higher redshifts 
to earlier times in the source frame, at which the afterglow is 


The progenitors of GRBs are believed to be stellar mass 
objects. So the intrinsic luminosity of GRB should not 
depend on the mass of their host galaxy, which has small 
mass at high redshifts; 

Long GRBs triggered by the death of massive stars, 
which have been shown to be associat ed with supernovae 
( Stanek et alJl2003 : Hiorth et~ai]l2003 ). provide a comple¬ 
mentary technique for measuring star formation rate. They 
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also represent a unique probe of the initial mass function 
and the star formation of massive stars at very high red- 
shifts. GRBs offer the exciting opportunity to detect the 
Population III (Pop III) stars; 

• GRB afterglows have smooth continuum spectra that al¬ 
low one to extract the IGM absorption features. The metal 
absorption lines in the GRB spectra make GRBs powerful 
sources to study the metal enrichment history; 

• The clean red damping wings of GRBs make them ideal 
tools to study the reionization of IGM and ISM properties 
of their hosts. 


First, GRBs can serves as the complementary tools to mea¬ 
sure dark energy and cosmic expansion. Type la super¬ 
novae (SNe la) are now treated as ideal “standard candles” 
for purposes of Hubble diagram. In 1998, two teams study¬ 
ing distant SNe la discovered independent evidence that the 


expansion of the Univ erse is speeding up (IRiess et alJ 11998 


Perlmutter et al. 1999), which is attributed to the mysterious 


component — dark energy. The accelerated expansion of 
the universe has also been confirmed by several independent 
observations including those of the cosmic microwave back¬ 
ground (CMB) (Spergel et al. 20031). th e baryonic acoustic os¬ 


cillations (BAO) (lEisens tein et al. 2005), X-ray gas mass frac¬ 


tion i n galaxy clusters ([A lien et al. 2004), and Hubble param¬ 


eters (Jimenez et al. 2003[ I Wang a nd Wangi l2014al) . The stan¬ 


dard ACDM model fits the observational data well, but other 
dark energy models could not be ruled out because of the pre¬ 
cision of current data. The latest observations confirm that 
about 70% of the energy density of the present Universe con¬ 
sists of dark energy. The direct evidence for the current ac¬ 
celeration of the universe is related to the observation of lu¬ 
minosity distances of high redshift SNe la. The SNe la can 
be observed when accreting white dwarf stars exceed the mass 
of the Chandrasekhar limit and explode. Thus they can be 
treated as an ideal standard candle. But hitherto the high¬ 


est redshift of SNe la is 1.914 (Jones et al. 2013). GRBs are 


shift, filling the gap between SNe la and CMB ( 

Dai et al. 2004: 

Ghirlanda et al. 2004tJ[ Liang and Zhang 

2005; 

Ghirlanda et al. 

2006a; Wang and Dai 2006a; Wang et al. 

2007t 

Schaefer 2007; 

Capozziello et al. 

2012; 

Wang 

2012). Similarly to SNe la, it 


has been proposed to use correlations of GRBs between var- 


emission to standardize GRB energetics (i,e., Ghirlanda et al. 

2004a 

Xu et al. 

2005| Liang and Zhang 2006; 

Firmani et al. 

2005; 

Amat|j200(: 

; Amati et al. 2008; Liang et al. 

2008, 2010). 


The high-redshift (z > 6) star formation history (SFH) is 
important in many fields in astrophysics. Direct star forma¬ 
tion rate (SFR) measurement beyond the reach of present in¬ 
struments, particularly at the low part of the galaxy lumi¬ 
nosity function. Long GRBs triggered by the death of mas¬ 
sive stars, provide a complementary technique for measuring 
the S FR ( Totanil 1 997k Wiiers et all 19981: iBromm and Loeb 
2002 ). Because the lifetimes of massive stars are short, the 
SFR can be treated as their death rate. Recent Swift ob- 


(Daigneetal. 2006; Le and Dernier 

2007; 

Yuksel and Kistler 

2007; Salvaterra and Chincarini 2007; Guetta and Piran 

2007; 

Kistler et al. 

2008; Campisi et al. 201(| Salvaterra et al. 

2012). 


The SFR revealed by the high-redshift GRBs seems to be much 
higher than that obtained from high-redshift galaxy surveys 
( Kistler et al. 2008, 200*1 Wang and Dai 2009). An enhance¬ 


ment about (1 + zy (5 ~ 0.5 - 1.5) in the observed rate of 


GRBs compared to SFR has been fou nd (IKistler et alJ 12008 


Robe rtson and Ellisl 1201 2: Wang 2013). In order to explain 


this discrepancy, many models have been proposed, includ - 
ing cosmic metallicity evolution (E l2008f, [Pin et ahl 12010 ), 
superconducting cosmic st rings ( Cheng et al. l2010h,~evolving 


star i nitial mass function ( Wang and Dai 201 la : Xu and Weil 


2008), evolution of the luminosity function break of GRBs 


Vir gili et alJ 1201 lb I Yu et all 120121) . But some other studies 
claimed that there is no discrepancy betw een GRB rate and SFR 
( Elliott et akll2012 : Hao and Yuanll2013 ). From host galaxy ob¬ 
servations, long GRBs prefer to form in a low-metallicity en¬ 
viron ment ( Le Floc’h et al. 20031: Stanek et al. 2006 ; Levesque 


2014), which is also required by theoretical prediction. The 


mass loss of stars through winds is proportional to the metal¬ 
licity, so low-metallicity can reduce the loss of angular momen¬ 
tum. Some studies have argued that GRB progenitors must have 
a low metallicity ( Woosle v and Heget 2006t Meszarosl 2006; 
Langer and Normanl2006l) . Observations also show differences 


in the population of GRB host galaxies compared to ex pecta¬ 


tions for an unbiased star-formation trace r (ITanvir et al 
Fruchter et al. 1120061 Svensson et al.ll2010h . 


pecta- 

120041 


The metal enrichment history has several important 


consequences for s tructure formation (iMadau et all 12001 


Karlsson et al.l 2012). An early phase of metal injection may 


qualitatively change the character of star formation, from a 
high-mass (Pop III) mode to a normal, low-mass dominated 
(Pop I/II) one, if the en richment exceeds a ‘critical metal¬ 
licity ’ ofZcrit ~ 10 4 Zr. ; ( Bromm et al. 2001 : Schneider et alJ 


20021 2006). This mode transition has crucial implica- 


(Bromm and Loeb 2006; Campisi et al. 20 

1; de Souza et al. 

2011), for the cosmic reionization (Cen 2003 

Wyithe and Loeb 


20031: Furlanetto and Loebl 2005), and for the chemical abun 


dance of low-met a llicity stars dQ ian and Wasserburgl 12001 
Frebel et al. 20071 2009; Jumlinsorj" 2010 ). Absorption 


lines on the spectra of bright background sources, such as 
GRBs or quasars, are main sources of information about 
the chemical properties of high-redshift Universe ( Ohl 12002; 
Furlanetto and Loeb 20031 Qppenheimer et al. 20091) , These 


lines are due to absorption by metals in low-ionization stages 
which arise in the higher column-density gas associated with 
damped Lya absorbers (DLAs). GRBs as bright sources have 
a number of a dvantages compared to traditional lighthouses 
such as quasars ( Bromm and Loebl2012 ). Their number density 
drops much less precipitously than quasars at z > 6 ( Fan et al.l 
l 2006l) . together with the power-law character of their spectra, 
renders them ideal probes of the early IGM. Quasars show 
strong spectral features, such as broad emission lines or the 
so-called “blue bump” that complicate the extraction of IGM 
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absorption features. The analysis of the spectrum of distant 
GRB 050904 has offered a wealth of detailed insight into the 
physical conditions of the host galaxy at z - 6.3 (ITotani et al. 
2006i) . Salvaterra et al. (2009) claimed that they identi¬ 


fied two absorption lines (Si IV and Fe II) in the spectrum 
of GRB 090423 , the most distant spectroscopically confi rmed 
burst at z = 8.2 (Salvaterra et al.112009 : Tanvir et al. 2009). 

Besides the above fields, high-redshift GRBs also may 
be useful tools to study dark matter and primordial non- 
Gaussianities. Because the matter power spectrum is dependent 
on the mass of dark matter, i.e., cold, warm and hot, which can 
result in different structure formation. So the number counts of 
high-redshift GR Bs can set strong lower limits on the dark mat¬ 
ter mass ( Mesinger et al. 20051: de Souza et al. 2013 ). The rate 
of high-redshift GRBs is also dependent on the amount of pri¬ 
mordial non-Gaussianity in the density field ( Maio et al]l2012h . 

In this work, we review the cosmological implications of 
GRBs. Then in the following sections the gamma-ray bursts 
cosmology is reviewed: The second section is dedicated to lu¬ 
minosity correlations and cosmological constraints from GRBs. 
Section 3 discusses the capability of GRBs to reveal the high- 
redshift SFR and cosmic reionization. In section 4, the capabil¬ 
ity of GRBs to probe the metal enrichment history is discussed. 
The last section provides a summary and future prospect. 


2. GRBs as standard candles to probe dark energy 

The best way to measure properties of the dark energy is to 
measure the expansion history of our universe, i.e., the redshift- 
distance relation. To this end, WFIRST has been proposed to 
determine the distances of 1000 SNe la with exquisite accuracy. 
To explore the properties of dark energy, the best method is to 
measure it over a wide range of redshifts, but SNe la can only 
be detected at low redshifts, i.e., z < 2.0. GRBs can extend the 
Hubble diagram to high redshifts. So many attempts have been 
performed to standardize GRBs. Frail et al. (2001) found that 
the collimated energetics of GRBs clustered around 5 x 10 50 


erg, which was confirmed by Bloom et al. ( 20031) . Observa¬ 


tions _aHojequiretihat_theGRB_ermssionJsonly in a small an¬ 
gle ( Waxman et al. 1998 : Fruchter et alJll999h . The collimated 
GRB model predicts that the a ppearance of an achroma tic break 
in the afterglow light curve ( iRhoads! 1997 ; Sarii [ 1999 ). This 


break is important to standardize the energetic of GRBs. 

Similar as SNe la, the luminosity correlations are required 
to probe dark energy using GRBs. In this section, we first re¬ 
view luminosity correlations of GRBs. Then the progress on 
dark energy revealed by GR Bs is discussed. Some reviews have 
discussed this topic (i.e. , Ghirlanda et alj|2006at Dai and W ang 


Lino tuple, yi.v-., v_ij.ij.jl lcijjviti v-c cu. v/u-, j-^cij tiuvi r * ci 

2007t Capozziello et al.luOl 2t Amati and Della Vallell20 13 ). 


2.1. The luminosity correlations of GRBs 

The luminosity correlations are connections between param¬ 
eters of the light curves and/or spectra with the GRB luminosity 
or energy. The isotropic luminosity can be calculated as 


l- — And, / J hoio 


( 1 ) 


and the total collimation-corrected energy is 

By — B'isoF'beam — And L S boio^ 1 b earn (1+Z). (2) 

Here, / J ho io and S boio are the bolometric peak flux and fluence, 
respectively, while Abeam = 1 — cos 0j et is the beaming factor 
with jet opening angle 0 ]ci . The peak fluxes and fluences are 
given over a wide variety of observed bandpasses, and with the 
wide range of redshifts which correspond to different range of 
energy bandsjn the rest frame of GRB. So the K-correction is 
important ( Bloom et al. 200 lh . / J ho |„ and .S’ boio are computed 
from the differential energy spectrum Off;) as follows: 


Rbolo 


fl 


10 4 /( 1 +:) 

1/(1+*) 


E<t>(E)dE 


f E Q>(E)dE 


> boio 


= S x 


r 10 4 /(J+z) 

Ji /(i +z) E®(E)dE 


J^ max E®(E)dE 


(3) 


(4) 


with P and S being the observed peak energy and fluence 
in units of photons/cm 2 /s and erg/cm 2 , respectively, and 
(Zimin, £max) the detection thresholds of the observing instru¬ 
ment. In general, the differential energy spectrum is modeled 
using a broken power - law ( Band et al.lll993l) . 


d>(£) = 


AE a e -(2 +a) E/E^ E < £^£ peak 


BE? 


(5) 


otherwise 


where a is the power-law index for photon energies below the 
break and (5 is the power-law index for photon energies above 
the break. Some differential energy spectra of GRBs also show 
power-law spectra plus an exponential cutoff. The luminosity 
distance di can be expressed as 


ddz) = 0 + z) 


- f 

H 0 Jo 


dz! 

E(z') 


where E 2 (z) = Qm( 1 +t:) 3 + G x / x (z) and fdz) is given by 


fdz) = exp 


1 f — 

Jo 1 


+ w(z) 


+ z 


d? 


( 6 ) 


(7) 


where w(z) is the equation of state (EOS) of dark energy. For 
ACDM, Eq.(|6| reduces to 


d L {z) = (\ + z) 


k£ 


dz' 


EIm( 1 + z0 3 ) + 


( 8 ) 


The promising luminosity correlations of GRBs are as fol¬ 
lows: 

• Lj so - riag correlation. The luminosity-time lag correla 


tion was first discovered by Norris et aD ( 2000h based on 
six GRBs observed by BATSE with optical redshifts (also 


see 


Schaefer et al. 


200 l j), which was confirmed by GRBs 


observed by Swift (Gehrels et al. 2006) including the pecu 


liar long GRB 060218 (Liang and/Zhang 


correlation is challenged by recent study (Bemardmi et al. 


2006). But this 
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Figure 1: The Amati correlation of 101 GRBs. Solid line is the best fit , and 
dashed lines represent the ltr dispersion. (Adapted from Figure 1 in lWang et alj 

1201 ilk ) 


2015). This correlation shows that more luminous bursts 

are also characterized by shorter time lags, i.e., Lj s 
—-1,25 

lag 


(Norris et al. 2000l) . This correlation has been used 


as a redshift indicator to estimate z for GRBs (Band et al 


2004), and to constrain cosmological pa rameters and dark 
energy ( SchaefeJl2007 : Wang et aPl2007 ). The interpreta¬ 
tion of this correlation is various, including viewing angle 
of the collimated i et dloka an d Nakamura 2001) or radia¬ 
tive cooling effect ( Crider et al. 119991) . 


Tiso ~ V correlation. Fenimore and Ramirez-Ruiz (2000) 
found that the time variability is correlated with the lumi¬ 
nosity of GRBs, which indicates that more luminous bursts 
have more variable light curve s. Later, this correlati on is 
confirmed with more GRBs bv iReichart et al.1 (|2001) and 
Schaefer et al. (2001 ). But the intrinsic scatter of Lj s n - V 
is very large ( Guidorzi et al. 2006 ; Wang_eUil. 201 lh. an d 
the index is still under debate did and Paczvnskil 20061) . 
The variability V is different for various instruments. The 
origin of the L ISI) - V correlation may be _based on the 
screening effect of the photosphere ( Kobavashi et al.l2002 ; 
Meszaros et al.ll2002l) . 


Amati correlation. Amati et al (2002) found that the 
isotropic energy Ej so is correlated with the rest-frame peak 
energy of the prompt spectrum, i.e., £p e ak K Sub- 


and Fermi ) confirmed this correlation ( 

Ama 

et al. 2009; 

Ghirlanda et al. 

2010; Sakamoto et al. 

2011; 

Wang et al. 


2011). Moreover, it was found that the Amati correlation 
also holds within individual GRBs using time-resolved 
spectra, and the slopes are consistent with the correla- 
tion from time-integ rated spectra ( Ghirlanda et al. 2010fc 
Frontera et al. 2012). The possible interpretations of Am¬ 
ati correlation include the synchrotron mechanism in rel¬ 
ativistic shocks (Zhang and Me szar os 2002) and emis¬ 


sion from off-axis relativist ic jets dYamazaki et al. 2004; 
Eich ler and Levinsonl 2004 ). Figure |T] shows the Amati 
correlation including 101 GRBs. The intrinsic scatter is 
0.62. 


Yonetoku correlation. The correlation A 


E 2 

- peak 


was found with a sample of 16 GRBs (Yoneto ku et al 
2004t Wei and G aol 2003 ). This correlation was con 


firmed by Liang et al. (2004). Similar as Amati correla¬ 
tion, £peak - Tiso also holds within individual GRBs using 
time-resolved spectra (Ghirl anda et al.ll2010l) . The possi¬ 
ble origin of this correlation is similar as that of Amati cor¬ 
relation. Figure [2] shows the Yonetoku correlation. The in¬ 
trinsic scatter of this correlation is 0.62 ( Wang et al . 2~0nj ). 


Ghirlanda correlation. A tight correlation between spec¬ 
tral peak ener gy and c o llimate d energy E y was 


discovered by Ghirlan da et al.l (I2004al) using 15 GRBs. 
The intrinsic scatter is up to 0.1, so this correlation is a 
promising tool to cons train dark energy (IDai et al.l 12004; 


Ghirlanda et al.l l2004bl). Bv consi dering the wind cir- 


cumburst density, Nava et al . (2006) found this correla 
tion also holds, and the intrinsic scatter is even smaller. 
One of the major challenge for this correlation is that 
most of GRBs observed by Swift do not show achro 


matic breaks in the afterglow light curve dWillingale et al 
2007). So the break time is very hard to determine. 


This correlation can be understood within the annular 
jet mo del dE ichler a nd Levinson 20061) and photosphere 
model (Thompson 2006|). The latest Ghirlanda correlation 
is shown in Figure [3] 


Lj so - £peak - 7().45 correlation. Firmani et al. ( 20061) found 
a tight correlation Lj so oc using parameters in 

prompt emission only. To .45 represents the “high signal” 
timescale. 


Liang-Zhang correlation. Without im posing any theo¬ 
retical model, Liang and Zhang (2005 0 found an empiri¬ 
cal correlation among the isotropic energy of the prompt 
gamma-ray emission Tj s0 , the rest-frame peak energy 
Epeak, the rest-frame break time in the optical band (break 
using 15 bursts. The correlation reads E pea k K E^o 2 /{^ k . 
If we take the optical break time as the jet break time, this 
correlation is similar to Ghirlanda correlation. The intrin¬ 
sic scatter of this correlation is als o very small, so it could 
be used to constrain dark energy ( Liang and Zhangl 2005 ; 
Wang and Dai 2006al; IWei et alil2013 ). The Liang-Zhang 
correlation is shown in Figure [4] 


Lx - T a correlation. iDainotti et al.1 (2 008 ) discovered a 
tight correlation between X-ray luminosity Lx and T a , 
where T a is the time at which the X-ray li ght c urve es¬ 


tablishes a afterglow power-law decay (IWillingale et al 
2007). The intrinsic scatter of this correlation is about 
0.33. By adding_a third parameter isotropic energy E VM , 


Xu and Huang (2012) found a new correlation, i.e.. Lx 

T’-0.87£0.88 
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Figure 2: The Yonetoku correlation. Solid line is the best fit, and dashed lines 
represent the lcr dispersion. (Adapted from Figure 1 in lWang et al.l d201 1.). t 



Figure 3: The Ghirlanda correlation. Solid line is the best fit, and dashed lines 
represent the lcr dispersion. (Adapted from Figure 1 in lWang et al.l d201 lid 



Figure 4: The Liang-Zhang correlation of 33 GRBs in the ACDM model. 
(Adapted from Figure 1 in Wei et al] 1201T ).) 


Due to GRBs cover large redshift range, whether the cor¬ 
relations evolve with the redshift should be discussed. It’s 
found that the slope of Amati correlation may vary with red¬ 
shift significantl y us ing a sm all sample of GRBs (Li 2007b . 
Basilakos and Perivolaropoulos (2008) found no statistically 
significant evidence for redshift dependence of slopes in five 
correlations using 69 GRBs compiled by Schaefer! ( 2007b . 
Wan g et al.l ( 2011 ) enlarge the GRB sample and test six GRB 


correlations. There is no statistically significant evidence for 
the evolution of the luminosity correlations with redshift is 
found. The slopes of correlations versus redshift are all con¬ 
sistent with zero at the 2 cr confidence level. 

The instrumental selection effects may affect the observed lu¬ 
minosity correlations. Some outliers to these correlations have 
been discovered (Soderberg et al. 2004; Vaughan et al. 2006; 
Campana et al. 2007; Rizzuto et al. 2007; Urata et al. 2009). 
Nakar and Piran (2005) considered the different samples of 
GRBs detected by BATSE. They found that a large fraction of 
GRBs were inconsistent with the Amati correlation by assum¬ 
ing all possible redshifts. Band and Preece (2005) also found 
that about 88 and 1.6 per cent of their BATSE GRBs were out¬ 
liers to the Amati and Ghirlanda correlations, respectively. But¬ 
ler et al. (2007) claimed that the Amati correlation exists, but 
it may be due to the selection effect. Using a large sample of 
GRBs with pseudo-redshifts determined by the L lso - T] ag corre¬ 
lation, Ghirlanda et al. (2005) argued that the Amati correlation 
really existed by taking into account the intrinsic scatter. By 
considering the triggering threshold limits for several GRB de¬ 
tectors, Nava et al. (2008) and Ghirlanda et al. (2008) claimed 
that only 6% of BATSE long GRBs are inconsistent with Amati 
correlation, while on outliers to the Ghirlanda correlation are 
found. By simulating the BATSE Large Area Detectors trig¬ 
gering thresholds, Shahmoradi and Nemiroff (2011) found that 
the Amati and Ghirlanda correlations are statistically real but 
strongly affected by the thresholds of GRB detectors. Ghirlanda 
et al. (2012) studied the selection effect on the Yonetoku cor¬ 
relation using Monte Carlo simulations, and found this correla¬ 
tion must be physical. Dainotti et al. (2014) proposed a general 
method to check the selection effects for GRB correlations and 
found the L x - T a correlation is not generated by the biases. In¬ 
terestingly, using the time-resolved spectra, similar correlations 
were found in individual bursts (Firmani et al. 2009; Ghirlanda 
et al. 2010; 2011). This strongly supports that the correlations 
are physical. 


2.2. Constraints on dark energy and cosmological parameters 

The most common method to constrain dark energy is 
through its influence on the expansion history of the universe, 
which can be extracted from the luminosity distance di (z) and 
the angular diameter distance <7,4(z). In addition, the weak grav¬ 
itational lensing, growth of large-scale structure, and redshift 
space distortion can also provide useful constraints on dark en¬ 
ergy. Theoretical models can be tested using the \ 2 statistic. 
The typical way to probe dark energy from standard candles 
is as follows. With luminosity distance di in units of mega- 
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parsecs, the theoretically predicted distance modulus is 
p = 5 logW/j + 25. 


(9) 


The likelihood functions for the cosmological parameters can 
be determined from statistics. 


X 


(Q. m ,Qde) = ^ —■—- 


lPi(Zi, Hq, D m , D 0£ ) - Ato,/]" 


( 10 ) 


W>j 


where /vo,/ is the observed distance modulus, and rr ;Jo . is the 
standard deviation in the individual distance modulus. The con¬ 
fidence regions in the Dm - Q.de plane can be found through 
marginalizing the likelihood functions over //q (i.e., integrating 
the probability density p <x exp(-y 2 /2) for all values of //(>). 

A lot of effort had been made to constrain cosmological 
parameters using GRBs_since their cosmological origin was 
confirmed. Schaefer ( 2003) obtained the first GRB Hubble dia¬ 
gram based on L ls() - V correlation, and found the mass density 
Qm < 0.35 at the lcr confidence level. After IGhirlanda et al 


(2004a) found the Ghirlanda correlation, Dai et al. (2004) 


first used this correlation with 12 bursts and found the mass 
density £1 m = 0.35 ± 0.15 at the lcr confident level for a flat 
universe by assuming that some physical explanation comes 


into existence. 


Ghirl anda e t al.J (2004b) using 14 GRBs and 
SNe la obtained Dm = 0.37 + 0.10 and Qa = 0.87 + 0.23. 
Assuming a flat universe, the cosmological parameters were 
constrained to be fly t — 0.29 ± 0.04 and Da = 0.71 ± 0.05 


(Ghirlanda et al. 2004b). For a flat universe, iFirmani et al 


( 20051) found Dm = 0.28 ± 0.03 and zt - 0.73 ± 0.09 for the 
combined GRB and SNe la sample, where zt is the transition 
redshift, at which the universe starts accelerating expansion. In 
the dark energy model of w z = wo, they found = 0.44 and 
wo = -1.68 with zt — 0-40 for the combined GRB and SNe la 


sample. Xu et ah (2005) obtained = 0.15*Jj'^(lc) using 17 


GRBs. Friedman & Bloom (2005) discussed several possible 
sources of systematic errors in using GRBs as standard candles. 
Using the /f pca k - £) so - /break correlation, Liang and Zhand 
( 2005 ) found the lcr constraints are 0.13 < Dm < 0.49 and 
0.50 < Da < 0.85 for a flat universe. They also found the 


Wang and Dai (2006a) using the 


transition redshift is 0.78/;g 22 . 

Liang-Zhang correlation to investigate the transition redshifts 
in different dark energy models via GRBs and SNe la, see also 
Wang and Dai ( 2006b ). Di Girolamo et al. (2005) simulated 
different samples of gamma-ray bursts and found that Dm could 
be determined with accuracy ~ 7% with data from 300 GRBs. 
Mea nwhile, many works ha ve been don e on this field, such 
as Mortsell and Sollermanl (2005), Bertolami and Silva 


(2006), 


Firmani et al 


(2006), Hooper and Dodelson 


< 2007 fj _ Lietakl ~ 2008al). iMosauera Cuesta et al 


BasilakosandPerirolaroEOulos d2008 ), _ Wang[ 


), Yuetal. (2009 

, Cardone et al. ( 

2009ft, 


Cardone et al. ( 

2011), 

Demianski et al. 


( 120131) . More recently, theoretical arguments and observational 


evidence both suggested that a small fraction of fast radio 
bursts (FRBs) may be associated with GRBs (lBannister et al 


2012i iThornton et al.1120 1 3t IZhangl 120141) . So the dispersion 
measure from FRBs and redshifts from GRBs makes these 
systems a plausible tool to study cosmological param eters 


( Deng and Zhangf2014 : Zhou et al. 2014 : Gao et al. 12014 ). 

Unfortunately, because of lack of low-z GRBs, the luminos¬ 
ity correlations has been obtained only from moderate-z GRBs. 
So this correlation is cosmology-dependent, i.e., the isotropic 
energy £) so , collimation-corrected energy E y , and luminosity 
Lj so are as functions of luminosity distance c// , which is de¬ 
pendent on cosmology model. This is the so-called “circularity 
problem” of GRBs. In the following, we present different meth¬ 
ods to overcome this problem. 

The first method is fitting the cosm ological parameters and 
luminosity correlation simultaneously. IGhirlanda et al.l (12004b) 
first used this method to overcome the “circularity problem”. 
Schaefer (2007) used 69 GRBs and five correlations to con¬ 
strain cosmological parameters. Wang et al. ( 20071 ) also used 
69 GRBs and other cosmological probes to constrain cosmo¬ 
logical parameters. They make simultaneous uses of five lumi¬ 
nosity indicators, which are correlations of L IS[) - ri ag , T ls „ - V, 
Bpeak - Liso, £peak - E y , and r RT - L iso . After obtaining the dis¬ 
tance modulus of each burst using one of these correlations, the 
real distance modulus can be calculated, 


A/fit 


= (Zw/^)/(ZO. 


( 11 ) 


where the summation runs from 1-5 over the correlations with 
available data, /r, is the best estimated distance modulus from 
the i-th relation, and <r lh is the corresponding uncertainty. The 
uncertainty of the distance modulus for each burst is 




(T M: y 


1/2 


( 12 ) 


Because each correlation with its limitations and possible bi¬ 
ases, simultaneous using may affect the cosmological indica¬ 
tors. But in order to enlarge the sample, some attempts have 
been performed. When calculating constraints on cosmologi¬ 
cal parameters and dark energy, the normalizations and slopes 
of the five correlations are marginalized. The marginalization 
method is to integrate over some parameters for all of its possi¬ 
ble values. The^ 2 value is 


2 

Tgrb 


[piiZi, Ho, D m , D 0£ ) - 


i=i 


cr„ 


(13) 


where p.(wj and cr mi are the fitted distance modulus and its error. 
In addition to GRBs, SNe la, CMB, BAO, X-ray gas mass frac¬ 
tion in galaxy clusters and growth rate data are also ideal cos¬ 
mological probes. CMB is the remnant of the cosmic recombi¬ 
nation epoch. It contains abundant information of the early uni¬ 
verse. The positions of the acoustic peaks contains the informa- 


Bond and Efstathiou 


cjpeaks i 

tionof dark energy dPeebles and Yulll970l : _ 

1984). The shift paramet er is defined as ( Bond et al. 1 19971 : 
Wan g and Mukherieell2006l) 


VDm . 

K = sinn 

vmn 


/|D,I I I = 1.70 ±0.03, (14) 

Jo E(z)J 
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where E(z) = H(z)/Hq and the function sinn(x) is defined as 
sinn(x) = sin(x) for a closed universe, sinn(x) = sinh(x) for 
an open universe and sinn(x) = x for a flat universe. The last 
scattering redshift zi s can be fitted as 

Zis = 1048 [ 1 + 0.00124(Qi/rr a738 ][l + gl (Q. M h 2 ) 82 ], (15) 


where the quantities gi and gz are defined as gi = 
0.078(Q fe /t 2 r a238 [l + 39.5(Q. h h 2 ) 0163 ]- 1 and g? = 0.5 6[1 + 
21.1(Qfe/i 2 ) 1-81 ] -1 respectively ( Hu and Sugivamaj 199 6). The 
X 2 value is 

2 CR - 1.70) 2 
Xcmb ~ o.03 2 


(16) 


BAO refers to an overdensity of baryonic matter due to 
acoustic waves which prop agated in the early universe (Silk 
1968; Peebles and Yuli 1970h . BAO provides a “standard ruler” 


for length scale in cosmology to explore the expansion his¬ 
tory of the universe. The acoustic signatures in the large-scale 
clustering of galaxies can be used to constrain cosmological 
param eters by detection of a peak in the correlation function 
( Eisenstein et al.ll2005l) . The A parameter is defined as 



A = 




Z 1 


Z 1 


1 


E(z i) |£2*| 


- simr 


dz 

Wz) 


1/3 


(17) 


measured from the SDSS data to be A =0.469(0.95/0.98) 0 35 + 
0.017, where zi = 0.35. The^- 2 value is 


Figure 5: Joint confidence intervals of lcr for (Qm, FI a) from the observational 
datasets. The thick black line contour, the blue contour, the red contour, the 
yellow contour, the violet contour, the orange contour, and the purple contour 
show constraints from from all the datasets, 26 galaxy clusters, 69 GRBs, CMB 
shift parameter, 182 SNe la, BAO, and 2dF Galaxy Redshift Survey, respec- 
tively. The t hin so lid line represents a flat universe. (Adapted from Figure 2 in 
Wang et all f2007l) .) 


Xbao ~ 


(A - 0.469) 2 
0.017 2 


08) 


In Figure [5] the constraint on the ACDM model is shown. 
Different color contours represent constraints from different 
data. The best fitted parameters are consistent with a flat ge¬ 
ometry. Figure [6] shows the cosmological constrain on the 
w = wo model from SNe la (blue) and GRBs (red). The 
combined constraint is shown as solid contours. The con¬ 


straint is much tighter by adding the GRB data. Fi et al 


( 2008b ) also performed a Markov Chain Monte Carlo (MCMC) 
global fitting analysis to overcome the circularity problem. The 
Ghirlanda correlation and 27 GRBs are used. They treated 
the slopes of Ghirlanda correlation and cosmological param¬ 
eters as free parameters and determine them simultaneously 
through MCMC analysis on GRB data together with other ob¬ 
servational data, such as SNe la, CMB and large-scale structure 
(FSS). Amati et al. ( 20081) measured the cosmological param¬ 
eters using Amati correlation using global fitting method. The 
extrinsic scatter was assumed on the parameter of // nC ak. but the 
cosmological-dependent value is Ei SO ( Ghirlanckfeoool) . 

The second method is to calibrate the correlations of GRBs 
using SNe la data at low redshifts. The principle of this method 
is that objects at the same redshift should have the same lumi¬ 
nosity distance in any cosmology model. Therefore, the lumi¬ 
nosity distance at any redshift in the redshift range of GRBs can 
be obtained by interpolating (or by other approaches) directly 
from the SNe la Hubble diagram. Then if further assuming 
these calibrated GRB correlations to be valid for all long GRB 




Figure 6: Joint confidence intervals of lcr to 3cr for {Q.m, wo) from the obser¬ 
vational datasets. The solid contours, the dashed blue contours, and the dot- 
dashed red contours show constraints from all the dat asets, 26 galaxy clu sters, 
and 69 GRBs, respectively. (Adapted from Figure 4 in lWang et alj fc007l) .) 
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data, the standard Hubble diagram method to constrain the cos¬ 
mological parameters from the GRB data_at high redshifts ob¬ 
tained by utilizing the correlations. Liang et al, (2008) first cali¬ 
brated the GRB correlations using an interpolation method. The 
error of distance modulus of linear in terpolation can be calcu¬ 
lated as ( Liang et aj]l2008 : Wei 11201 (ih 


Vn = ([(z;+i - z)/(z i+ i - Zi )] 2 e£; + t(z - Zi)/(z>+i - z,)] 2 f^ ;+1 ) 


1/2 


where e fl j and e, 


(19) 

-,/jj+i are errors of the SNe la, fJ, and p,, | are the 
distance moduli of the SNe la at z, and z,+ 1 re specti vely. Similar 
to the interpolation method, Cardo ne et al. (2009) constructed 
an updated GRBs Hubble diagram on six correlations calibrated 
by local regression from SNe la. Kodama et al. (2008) pre¬ 
sented that the Li so - £,’ pea k correlation can be calibrated with the 
empirical formula fitted from the luminosity distance of SNe la. 
This method has been used to constrain cosmological parame¬ 
ters by c omb ining these GRB data with SNe la in a following 


work by Tsutsui et al. (2008). 


However, it must be noted that this calibration procedure de¬ 
pends seriously on the choice of the formula and various possi¬ 
ble formulas can be fitted from the SNe la data that could give 
different calibration results of GRBs. As the cosmological con¬ 
straints from GRBs are sensitive to GRBs calibration results 
( Wang !2008b . the reliability of this method should be tested 


carefully. Moreover, as pointed out by Wang (2008), the GRB 


luminosity correlations which are calibrated by this way are no 
longer completely independent of all the SNe la data points. 
Therefore these GRB data can not be used to directly combine 
with the whole SNe la dataset to constrain cosmological param- 


the fitting formula, 

Wang et al. ( 

2009a) used the cosmographic 

parameters (Canozzicllo and Izzo 2008; Vitagliano et al. 2010c 

Xia et al. 2012; Gao et al. Ebll 

. The luminosity distance can 

be expanded as (Visse^2004 

) 


1 a 

Cl ~ ~7Fa’ 

(21) 

1 Cl 

7 = IF a’ 

(22) 

1 a 

S ~ IF a' 

(23) 


Equation (l20t is only dependent on the cosmological principles 
and FRW metric, so the expansion is model-independent. But 
the Taylor-expansion of d/ is not valid at z > L So expansion d L 
as a function of y = z! 1 + z is much useful ( Cattoen and VisseJ 
2007t Wang and Dail20lTbh . The calibrated Hubble diagram of 
GRBs is shown in Figure I71( Wan g and DafeOl lbh . 



Figure 7: The Hubble diagram of 557 S Ne la (blue) and 66 high-redshift GRBs 
(red). (Adapted from Figure 2 in lWang and Dai (2011b ).) 


The third method is to calibrate the standard candles us 
ing GR Bs in a narrow redshift range ( 6z ) near a fiducial red 


shift JLa mb et al J2005 : Liang and Zhand2006l : Ghirlanda et al.. 

2006b). Liang and Zhang ( 20061) proposed a procedure to cali¬ 
brate the Liang-Zhang correlation with a sample of GRBs in a 
narrow redshift range. No low-redshift GRB sample is needed 
in this method. The calibration procedure can be described as 
follows. First, calibrate the power-law index of Liang-Zhang 
correlation using a sample of GRBs that satisfy this correlation 
and are distributed in a narrow redshift range. The power-law 
index can be derived using a multiple regression method. Sec¬ 
ond, marginalize the coefficient value over a reasonable range. 

However, the gravitational lensing by random fluctuations in 
the interve ning matter distribution induces a dispersio n in GRB 


dL = 7^{ Z + ~ qo)zl ~ \ ( : ~~ qo ~ + j0 ) z3 

+ [2 - 2q 0 - 15<75 “ 15 < 7 o + 5/o + 10q 0 jo + s 0 ] z 4 + 0(z 5 )}, 

where q is the deceleration parameter, j is the so-called “jerk”, 
and s is the so-called “snap” parameter. These quantities are 
defined as 


the intervening matter distribution induces a dispersion in ORB 
brightness ( Qguri and Takahashil2006l : Schaefej2007 ). degrad¬ 


ing their value as standard candles as well as SNe la (Hoi: 


119981) . GRBs can be magnified (or reduced) by the gravitational 
tensing produced by the structure of the Universe. The gravi¬ 
tational lensing has sometimes a great impact on high-redshift 
GRBs. First, the probability distribution functions (PDFs) of 
gravitational lensing magnification have much higher disper- 
sions and are markedly different from the Gaussian distribu¬ 


tion (ValageasIbOOPtlOguri and Takahashil20()5 Wang and Pail 

201 lbh . Figure l8l shows the magnification probability distri¬ 
bution functions of gravitational lensing at different redshifts 
( Wang and Daill201 lbh . Second, there is effectively a threshold 
for the detection in the burst apparent brightness. With gravi¬ 
tational lensing, bursts just below this threshold might be mag¬ 
nified in brightness and detected, whereas bursts just beyond 
this threshold might be reduced in brightness and excluded. 
Wang and Dai (2011b) considered the weak lensing effect on 
cosmological parameters derived from GRBs, and found that 
the most probable value of the observed matter density Om is 
slightly lower than its actual value, see Figure [9] The weak 
gravitational lensing also affects the dark energy equation of 
state by shifting it to a more negative value. 
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Figure 8: Magnification probability distribution functions of gra vitational lens- 
ing at re dshifts z — 1 , z = 3 and z = 7. (Adapted from Figure 5 in lWang and Pail 

feoilbl) .) 



redshift z 


Figure 10: Estimation of the uncorrelated dark energy EOS parameters at dif¬ 
ferent redshift bins (wi, W 2 ,tvio) from SNe Ia+BAO+WMAP9+H(z)+GRB 
data. The open points show the best fit value. The error bars are ltr and 2<r 
confidence levels. The dotted line shows the cosmological constant. (Adapted 
from Figure 3 in lWang and Pail 12014J) .) 


2.3. The equation of state of dark energy 



Figure 9: Confidence contours of likelihood from lcr to 3cr in the ACDM 
model. The black line contours from 116 GRBs and the dotted contours 
from 116 GRBs including magnification bias. (Adapted from Figure 7 in 
Wang and Dal j20lfbh .l 


The dark energy equation of state w is the most important 
parameter that describes the properties of dark energy. Whether 
and how it evolves with time is crucial for revealing the physics 
of dark energy. GRBs can provide the high-redshift evolution 
property of dark energy. The procedure is to bin w in z, and fit 
the w in each bin to observational data by assuming that w is 
constant in each bin. The function f(z) in equation (JTJ should 
be described as 


f(Zn -1 <Z<Zn) = ( 1 + Z) 3(l+ "'" ) f^ 1 + *> 

/=0 


3(Wj-w i+l ) 


(24) 


where w,- is the EOS parameter in the 1 th redshift bin defined by 
an upper boundary at Zi, and the zeroth bin is defined as zo = 0. 
Oi et al . (2008a) used GRBs and other cosmological observa¬ 
tions to construct evolution of the equation of state, and found 
that the equation of state w is c onsistent with the cosmological 
constant (also see Oi et al.l 2008b ). The confidence interval of 
the uncorrelated equation of state parameter can be significantly 
reduced by adding GRBs. After calib rating the GRB correla¬ 
tions using cosmographic parameters, IWang and Pail (1201 lbl) 
found that the high-redshift (1.4 < z < 8.2) equation of state 
is consistent with the cosmological constant. But some studies 
found that the equation of state w may deviate from -1 (i.e., 
Oi et al] 2009: ZhaoeUil. 2012h . In light of the Planck CMB 
data, Wang and Pail d2014ah found that the EOS is consistent 
with the cosmological constant at the 2 cr confidence level, not 
preferring to a dynamical dark energy, which is shown in Figure 

To] 
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Figure 11: Distribution of the isotropic-equivalent luminosity for 157 long- 
duration Swift GRBs. The shaded are a appro ximates the detection threshold of 
Swift BAT. (Adapted from Figure 1 in Wang 1 201 3l >. ) 


3. Probing the star formation rate 


3.1. Star formation rate derived from GRBs 

The association of long GRBs with core-collapse super- 
novae has been confirmed from obse rvations in recent years 
( Stanek et al. 2003 : Hiorth et~aj]l2003 ). which provides a com- 
plemen t ary te c hnique for measur i ng the high-redshif t SFR 
(Totani j_997^_ Wi jersjdal] 19981: Lamb a nd Reichartl 2000; 


Porciani and Madau 200 it iBromm and Loeb 2002 ). The se 


lection effects should be considered (for a review, see Coward 


2007). But one crucial problem appears, i.e., how to cali- 


may 

play an important role (Nataraian et al. 2005 

Daigne et al. 

2006 

Salvaterra and Chinearini 2007; Salvaterra et al. 2009; 

Camnisi et al. 

2010 ; 

Wanderman and Piran 2010; Cao et al. 

201JJ. Before the launch of Swift ( 

Gehrels et al. 

2004), the lu- 


minosity function is determined by fitting the observed log N 


log P distribution dSchmidl 1999c Porciani and Madaul 2001 


Guetta et alJl2005t iNataraian et al.ll2005l) . Thanks to the Swift, 
more redshifts of GRBs are measured. This makes it pos- 
sible to give more information on the luminosity functio n 
( Wanderman and Piran 201(lt Cao et al. 2011 : Tan et al. 2013 ). 


Because the form of luminosity function should be assumed 
and the model parameters of luminosity function is degener¬ 
ate, it is not easy to determine the luminosity function. A 
straightforward way to estimate the luminosity function is 


proposed by l.vnden -Bc 


_ _ - _ (119771) and then further developed 

by Efron and PetrosianTi 19921). This method has been used 
for GRB s ( Llvd-Ronning et all 2002 : Yonetoku et al. 2004 


Wu et al. 2012h . There are two luminosity function models in 
the literature, a broken power law and a single power law with 
an exponential cut-off at low luminosities. 

In order to avoid the poorly known luminosity function when 
studying high-redshift SFR, a method that only high-luminosity 
GRBs are used is proposed (lYiiksel et al. 2008: iKistler et al 


20091 1 Wang and Pail 120091: lYu et aklEoIi IWanglEoIl . The 
expected redshift distribution of GRBs is 


dN .s(z)p„(z) dV com /dz 
~r = F ^~F7 - 7 -;-’ 

dz (/beam) Z 


(25) 


where F(z) represents the ability to obtain the redshift, e(z) ac¬ 
counts for the fraction of stars producing GRBs, and p,(z.) is the 
SFR density. The F(z) can be treated as constant when we con¬ 
sider the bright bursts with luminosities sufficient to be detected 
within an entire redshift range. GRBs that are unobservable due 
to beaming are accounted for through (/beam)- The s(z) can be 
parameterized as e(z) = £o( 1 + z) s , where «o is an unknown 
constant that includes the absolute conversion from the SFR to 


the GRB rate in a given GRB luminosity range. Kistler et al 


(2008) found the index 6 = 1.5 from 63 Swift GRBs. A little 
smaller value 6 ~ 0.5 -1.2 has been inf erred from update Swift 
GRBs ( Kistler et al. 2009t Wang 20131) . In a flat universe, the 
comoving volume is calculated by 


dV com , r . o dD com 
- = D- om - 


dz 


dz 


where the comoving distance is 


B^com(z) 




dz' 


sjQni 1 + Z') 3 + 


(26) 


(27) 


In the calculations, the ACDM model with O,,, = 0.27, 
Q a = 0.73 and Hq-1\ km s 1 Mpc 1 from the Wilkinson Mi¬ 
crowave^ Anisotropy Probe (WMAP) seven-year data is used 
dKomatsu et al.lEoiH) . 

Figure [TT1 shows the isotropic luminosity distribution of 157 
Swift GRBs. The isotropic luminosity can be obtained by 


Tjso - Tj so ( 1 + z)/T 90 , 


(28) 


where Tgo is the duration time. The shaded area approximates 
the detection threshold of Swift BAT, which has a flux limit 
~ F\i m = 1.2 x 10~ 8 erg cm -2 s~*. So the selection effect is 
important. In order to exclude faint low -redshift GRBs that 
could not be visible at high redshifts, we only select luminous 
bursts. The luminosity cut L lv , > 10 51 erg s _1 is chosen in the 
redshift bin 0-4 (Yiiksel et al. 2008), which removes many 


low-redshift, low-luminosity bursts that could not be detected 
at higher redshift. The cumulative distribution of GRB redshift 
can be expressed as 


N(<z) _ N(0, z) 

Y(< Zmax) N(0, Ziriax ) 


(29) 


The value of z max is taken as 4.0. Because the SFR has been well 
measured at z < 4.0 ( Hopkins and Beacom 20061) . The theory 
predicted and observed cumulative GRB distributions is shown 
in Figure[12] The Kolmogorov-Smirnov statistic gives the min¬ 
imization for 6 = 0.5 ( Wangll2()T3 ). At the 2cr confidence level, 
the value of 6 is in the range -0.15 < 6 < 1.6. 

There are four redshift bins, z = 4 - 5, 5-6, 6-7, 7-8.5 and 
8.5 - 10. The GRBs in z = 1-4 play as a “control group” to 
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Figure 12: Cumulative distribution of 92 Swift GRBs with Lj s0 > 10 51 erg s -1 
in z = 0 - 4 (stepwise solid line). The dashed line shows the GRB rate inferred 
from the star formation history of Hopkins & Beacom (2006). The solid line 
shows the GRB rate inferred from the star fo rmation history including (1 + z) 0 5 
evolution. (Adapted from Figure 2 in ]Wand f2013l) .) 


Figure 13: The cosmic star formation history. The grey points are taken from 
Hopkins & Beacom (2006), the dashed line shows their fitting result. The trian¬ 
gular points_are from Bou wens e t al. (2009, 2011). The open circles are taken 
from iRob ertson a nd Ellisl i2012l) . The filled circles are the SFR derived from 
GRBs in Wang (2013). (Adapted from Figure 3 in lWand d2013l) .) 


constrain the GRB-to-SFR conversion. The theoretically predi¬ 
cated number of GRBs in this bin can be calculated as 


N\% = At 


= A 


A Q f 
4 n J| 

I 


dz F(z) s(z) 


P*(z) dV com /dz 

(/beam) 1 + Z 


dz( 1 +z) p*(z) 


s . / dVcom/dz 


1 + z 


(30) 


where A = At API Fo/47t(fi jeam ) depends on the total observed 
time of Swift, At, and the angular sky coverage, AO. The theo¬ 
retical number of GRBs in redshift bin zi - z .2 is 


N" 


/• \ a f 2 j 1 1 , \d dVcom/dz 
(P*) ZI -Z 2 A dz.(] + z) , (31) 

Jz i i Z 


where (p*) zl _ Z2 is the SFR in the redshift range zi - Zi . Repre¬ 
senting the predicated numbers, A*_ Z2 with the observed GRB 
counts, N z p 2 , we obtain the SFR in the redshift range zi - Zi, 




N^ Z2 f l 4 dz d -^(l + z) 6 Mz) 


N\ 


obs 

1-4 


pdz 


dVcom/dz 

1 +z 


(1 +z) 5 


(32) 


The derived SFR from GRBs are shown as filled circles in Fig¬ 
ure [13] Error bars correspond to 68% Poisson confidence in¬ 
tervals for the binned events ( Gehrelsl[l986 ). The high-redshift 
SFRs obviously decrease with increasing redshifts, although an 
oscillation may exist. The SFRs from GRBs are dramatically 
larger than those from other observations. The main reason is 
that other observations probe only the brightest galaxies, espe¬ 
cially at high redshifts. But GRBs can reveal the faint galax¬ 
ies at high redshifts due to their high luminosity. The SFR at 
Z > 4.48 is proportional to (1 + z) -3 , which is shown as solid 
line in Figure [13] 


3.2. Possible origins of high-redshift GRB rate excess 

Recent studies show that the rate of GRBs does not strictly 


high redshift (Le and Dermer 

2007; 

Salvaterra and Chine arini 

2007; Kistler et al. 2008; Yiiksel et al 

2008: Wang etal.l2009b: 

Robertson and Ellis 2012j Wan <1 

2013 

). The SFR inferred from 


the high-redshift (z > 6) GRBs seems to be too high in com¬ 
parison with the SFR obtained from some high-redshift galaxy 
i ( Bouwens et al.ll2009[ 2011 ). 


surveys 


3.2.1. Metallicity evolution 

A natural origin of the high-redshift GRB rate excess is the 
metallcity evolution. Theory and observation both support that 
long GRBs prefer to occurring in low-metallicity environment. 
Some theoretical studies of long GRBs progenitors using stel¬ 
lar evolution models suggest that low metallicity may be a nec¬ 
essary condition for a long GRB to occur. For popular col¬ 
lapse models of long GRBs, stars with masses > 3 0M 0 can 


be able to create a black hole (BH) remnant ( Wooslevl 1993 


Hirschi et al. 2005!) . The preservation of high angular momen¬ 
tum and high-stellar mass at the time of collapse ( Woosle v 
1993 ; MacFadven and Wooslev 20061) is crucial for producing 
a relativistic jet and high luminosity. Low-metallicity (0.1 - 
O.3Z 0 ) progenitors can theoretically retain more of their mass 
due to smaller line-driven ste llar winds (Ku dritzki and Puls 
2000t Vink and de Koterl 20051). and hence i preserve their an 


pr 

gular momentum (Yoon and Langer 2005; Yoon et al. 2006 


Woo slev and Hegei 20061) . because the wind-driven mass loss 


of massive stars is proportional to the metallicity. Observa¬ 
tions of long GRB host galaxies also show that they are typ¬ 
ically in low metal licity environment, for several local lon g 
GRB host galaxies ( Sollerman et al j 2006 : Stanek et al. 2006), 
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as well as in distant lo ng GRB hosts (i.e., Fruc hter et al.l l2006: 
IProchaska et al.lE007 ). 

□ d2008h studied the possibility of interpreting the ob¬ 
served discrepancy between the GRB rate history and the 
star formation rate history using cosmic metallicity evolution 
(Rustler et al. 2008). Under the assumption that the formation 
of long GRBs follows the cosmic star formation history and 
form preferentially in low-metallicity galaxies, the rate of GRB 
is given by 


Rgrb(z) = kGRB^lZth ,z)p*(z). 


(33) 


where &grb is the GRB formation efficiency, £(Z t h,z) is the 
fraction of galaxies at red shift z with metallicity below Z th 
( Langer and No rman 20 061) and p*(z) is the obs erved SFR. The 
function 2(Z t h,z) is ( Langer and Normanl2006l) 


UZ<h,z) = 


f[a 1+ 2,(Z th /Z Q ) 2 10 01 ^] 

T(ai + 2) 


(34) 



Figure 15: GRB host galaxy mass distribution. The solid lines represent the 
upper limits of the stellar mass of a GRB host galaxy given a metallicity cutoff 
of 12 + log(0 /H)kk 04 = 8.7 (black), and 12 + log(O/H)KK 04 = 8.6 (red). The 
dashed lines represent the lcr scatter. (Adapted from Figure 4 in lWang and Pail 
<20145) ,) 


where f and T are the incomplete and complete gamma func 


tions, a i 


= —1,16 and /3 = 2 ( Savaglio et al. 2005 ). IP 


(2008) found that the distribution of luminosity and cumu¬ 
lative distribution of redshift could be well fitted if Z± = 
0.3Z Q is adopted. Wang and Dai ( 20091) studied the high- 


redshift SFR by considering the GRBs tracing the star forma¬ 
tion history and the cosmic metallicity evolution. They found 
the SFR derived from GRBs is m arginal consistent with that 
from traditional way (i.e., Hopkins and Beacom 20061) . Us¬ 
ing Monte Carlo simulations. Pin et al.f ~( 201C ) compared the 
simulation results to the Swift observations with log N - log P 
and luminosity-redshift distributions. They found that the ob¬ 
served distributions are well consistent with that from simu¬ 
lations if the GRB rate is proportional to the SFR incorporat¬ 
ing with the cosmic metallicity history with Z,h = 0.6Z o . Fig¬ 
ure M shows the comparison between simulations and obser¬ 


vation. 


Wei et alJ ( 2014 ) examined the influence on the GRB 
distribution due to the background cosmology, i.e., R/, - ct 
Universe. However, a few GRB hosts with high metallicity 
are observed (i.e. GRB 020819), so that the role of metal- 


it is 

still debated (Price et al. 2007; Wolf and Podsiadlowski 

2007 

; Kocevski et al. 2009; Graham et al. 2009; Sve 

nsson et al. 

2010 

). For excellent reviews, see Fynbo et al. ( 

2012h and 

Levesque (2014). But there are some uncertainties when 


measure the metallicities of GRBs’ explosion region at_high- 
redsh i fts, such as chemical inhomogene ity dLevesaue et al. 
2010| Niino 2011 ). Wang and Dai ( 201 4bl) studied the metal¬ 


licity role from two aspects, the GRB host galaxies and red¬ 
shift distribution. They found that the the observed GRB 
host galaxy masses and the cumulative redshift distribution 
can fit the predicted distributions well if GRBs occur in low- 
metallicity 12 + log(O/H)KK 04 < 8.7, which is shown in Figure 
[T5l Trenti et al. (2015) found th at there is clear evidence for a 
relation between SFR and GRB ( Jimenez and Piranll2013l) . But 
a sharp cut-off of metallicity is ruled out. 


3.2.2. Evolving star initial mass function 

Wang and Pail ( 2011a ) proposed that the GRB rate excess 


may be due to the evolution of star initial mass function (IMF), 
also see (Xu and Wei 2008). Because an “top-heavy” IMF will 
lead to more massive stars at high-redshift which can result in 
much more GRBs. Considering long GRBs trace SFR, the rate 
of GRBs in an evolving IMF is 


Rgrb ' 


Nm>30M o 

V 


= K 


fin 


' 3 


0/ £" mf(m)d log / 


P*(z), (35) 


where K is a constant to be constrained and Rgrb is the rate 
of GRBs, representing the number of GRBs per unit ti me per 
unit volume at redshift z. The evolving IMF proposed by Dave 
d2008h is 


dN 

d log m 


■ oc 


m 0,3 for m < niiMF 

W _1 3 for III > 777lMF> 


(36) 


where ;«imf = 0.5(1 + z) 2 M Q , which has been constrained by 
requiring non-evolving star formation activity parameter. Fig- 
ure[T6]shows that the observed cumulative distribution of GRBs 
can be well produced by this model. 

3.2.3. Evolving luminosity function break 

Virgili et~afl ( 201 1) found that if the break of luminosity 
function evolves with redshift, the distributions of luminos¬ 
ity, redshift and peak photon flux from the BATSE and Swift 
data can be reproduced from simulations. The break luminosity 
function evolution can be in a moderate way oc L fc x(l +z)~ a8 ~'- 2 . 
Campisi et al. ( 2010 ) studied the luminosity function, the rate 
of long GRBs at high redshift, using high-resolution N-body 
simulations. A strongly evolving luminosity function with no 
metallicity cut may well explain the log N - log P distribution 
of BATSE and Swift data. 
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Figure 14: Comparisons of log L - logs distributions and log L, logs, and log P distributions between the observed Swift/BfiiT GRB sample (solid) and simulations 
(open dots and dashed lines) for: /?grb(z) k P*(z) x 2.(Z t h,z). left four panels are for the trigged GRB sample. Right four panels are for the sample with redshift 
measurement. One dimensional K-S test probabilities for the comparisons are presented in each panel. (Adapted from Figure 5 in lOin et all 120101) .) 



Figure 16: The cumulative distribution of 72 Swift long GRBs with /. lsn > 
0.8 X 10 SI erg s 1 (stepwise solid line). The dotted line shows the GRB rate 
inferred from the star formation history of Hopkins & Beacom (2006). The 
dashed line shows the GRB rate inferred from star formation history including 
an evolving IMF. (Adapted from Figure 2in lWang and Pail 1201 lall .) 


3.2.4. Superconducting cosmic string 

Cosmic strings are thought to be linear topological defects 
that could be formed at a phase transition in very early Uni¬ 
verse. By considering that high-redshift GRBs 080913 and 
090423 are electrom agneti c bursts of superconducting cosmic 
strings, Cheng et al. ( 2010h s howed the high-red shift GRB ex¬ 
cess can be reconciled. But Wang et al. ( 201 lh claimed that 
GRBs from cosmic string have a very small angle, about 10 ~ 3 , 
which could be in contradiction with the opening angle of the 
GRB outflow. Cheng et~ai] ( 2011 ) pointed out that the angle 
is not the opening angle of the GRB outflow, but is just the 
collimation angle of the radiation of the corresponding string 
segment. We must caution that the existence of cosmic string is 
only speculative. 


4. Probing the Pop III stars and High-Redshift IGM 

4.1. Observational signature of Pop Ill GRBs 

The first stars, also called Population III (Pop III) stars, are 
predicted to have formed in minihaloes with virial tempera¬ 
tures r vir < 10 4 K at z > 15 ( Tegmark et al.lll99~7 ; Yoshida et al 


20031 Bromm and Larson 2004 ). Numerical simulations show 
that Pop III stars forming in primordial minihaloes, were pre- 
domina ntly very massive stars with typical masses M t > 
1 00 M r, ( Bromm et al.l 1 99St Bromm and Loebl2002 ; Abeletal 


2002), for recent reviews, see Bromm et al. (2002) and Bromm 


( 2013 ). They had likely played a crucial role in early universe 
evolution, including reionization, metal enrichment history. 


called Pop III GRBs (Hegeretal 

2003; 

Bromm and Loeb 

2006; Komissarov and Barkov 

2010 

Stacv et al. 2011), which 


(Tomaetal. 2011; 

'lagakura et al. 

2011 ; 

Camnisi et al. 2011; 

Meszaros and Rees 

2010^ Nakauchi et al. 

2012h. Direct obser- 


vations of the Pop III stars have so far been out of reach. The 
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Figure 17: Left panel Minihalo circumburst density. Shown is the hydrogen number density as a function of distance from the central Pop III star at the moment 
of its death. Typical circumburst densities are ~ 1 cm -3 . Right panel Atomic-cooling-halo circumburst density. The density profiles are calculated from the Shu 
solution. T he cas e of photoheating from only a single Pop III star (solid lines), and that from a stellar cluster {dotted lines). (Adapted from Figures 1 and 2 in 
IWang et ni 5oii .) 


properties of Pop III stars may be revealed by their remanents. 
Pop III GRBs. 

In order to predict the observational signature of Pop III 
GRBs, the unusual circumburst environment that hosted the 
Pop III stars should be determined. In particular, the proper¬ 
ties of the afterglow emission of Pop III GRBs depend on the 
circumburst dens ity ( Ciardi and Loeb 2000l Gou et al. 2004 : 
Wan g et al. 2012h . In particular, the central minihalo envi¬ 
ronments just before a massive star die and a GRB bursts 
out can be understood as follows. The number of ionizing 
photons depends strongly on the central stellar mass, which 
is de termined by a accretion flow onto the growing protostar 
(e.g., McKee and Tanl 20081 Hosokawa et al. 2011 : Stacv et al 


2012). Meanwhile the accretion is also affected by this radia¬ 
tion field. So the assembly of the Pop III stars and the devel¬ 
opment of an H II region around them proceed simultaneously, 
and affect each other. The shallow potential wells of miniha¬ 
los are unable to maintain photo-ionized gas, so that the gas 
is effectively blown out of the minihalo. The resul ting photo¬ 
evaporation has been s tudied ( Alvarez et al. 20061: Abel et al 
2007tlGreifetai1l2009l) . 


The photoevaporation from minihalos can be described as the 
self-similar solution for a champagne flow (Shu et al. 2002). 
Assuming a p oc r~ 2 density profile, the spherically symmet¬ 
ric continuity and Euler equations for isothermal gas can be 
described as follows: 


[(v-x) 2 - 1] 


1 da 
a dx 


dv 


a -(x - v') 


(x - V'), 


[(v-x)" - 1]— = 
dx 


(x - v)or ■ 


(X-V), 


(37) 


(38) 


where x = r/c s t, and p{r,t) = a(x)/4nGt 2 = mun(r)/X and 
u(r, t) = Cjv(x) are the reduced density and velocity, respec¬ 
tively. c s is the sound speed and X = 0.75 the hydrogen mass 


fraction. We set the typical lifetime of a massive Pop III star as 
t = u - 3 x 10 6 yr. 

In the left panel of Figure [IT] we show the density profiles 
at the end of the Pop III progenitor’s life in the minihalo case. 
The circumstellar densities are nearly uniform at small radii. 
Such a flat density profile is markedly different from that cre¬ 
ated by stellar winds. But in the atomic cooling halo case, 
star formation and radiative feedback is not wel l understood 
( Johnson et al.ll2009t Safranek-Shrader e t all201 2h. such a s the 


masses of stars, and stellar multiplicity (Clark et al. 2011). So 
we also use the formalism of the Shu solution as above. We as¬ 
sume that either one Pop III star or a small stellar cluster forms. 
The densities are shown in the right panel of Figure [17] Similar 
to the minihalo case, densities are nearly constant at small radii, 
but overall values are much higher, which is due to the deep 
potential wells, so that photoheated gas can easily be retained. 
Typical circumburst densities are n ~ 100 cm -3 . Pop III GRBs 
originating in atomic cooling halos may be extremely bright. 

The 


typical parameters of the afterglow emission are 
adopted, T 0 = 300, E mo = 10 53 erg, Ao = 10 12 cm, e e = 0.3, 
£b = 0.1, and p = 2.5. As an example, in FigureQjO the M-band 
(v = 6.3 x 10 I3 Hz) light curve is shown. Figure [19] gives the 
observed flux at v = 1.36 x 10 14 Hz as a function of redshift 
in the minihalo case. The lines with filled dots, black trian¬ 
gles and open dots correspond to an observed time of 6 min¬ 
utes, 1 hour, and 1 day respectively. The straight line marks 
the K-band sensitivity for the near-infrare d spectrograph (NIR - 


Spec) on James Webb Space Telescope (Gardner et al. 2006!) . 


The high-redshift cut-off is due to the Lya absorption. The flux 
will be completely absorbed by the intervening neutral IGM. At 
these frequencies, the flux of afterglow is weakly dependent on 
redshift of GRB. There are two reasons. First, the time dilation 
effect implies that the high redshift means the earlier emission 
times , where the afterglow are much brighter ( Ciardi and Loebl 
2000; Bromm and Loeb 2012h . Second, circumburst densities 
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Figure 18: Light curve at v = 6.3 X 10 13 Hz (M band) of Pop III GRBs. The 
emission from the forward shock (dashed line), the reverse shock (solid red 
line), and their combination (solid black line) are shown. (Adapted from Figure 
5 in lWang et alj 120121) .) 


of GRBs modestly increase with redshift. 

4.2. Metal enrichment history 

The metal enrichment history has several important influ¬ 
ences for cosmic structure formation. For example, the metal 


inject i on change the mode of star formation (IBromm et al 


2001 : Schneider et al. 20021) . The transition between Pop III 
star formation and “normal” (Pop I/II) star formation has im¬ 
porta nt implications, e.g., the expected GRB redshift distribu- 
tion ( Bromm and Loebl 2002h . reionization ( Wvithe and Loebl 
20031) , and the chemical abundance patterns of stars. So it is im¬ 
portant to map the topology of pre-galactic metal enrichment. 
Ten or thirty meter-class telescopes have been proposed to mea¬ 
sure the z > 6 IGM metallicity with the GRB afterglow (Oh 
20021) . Meanwhile, the relative gas column density from metal 
absorption lines can refl ect the enrichment history (IHartmann 
2008tlWang etakll2012l) . 

Absorption processes and absorption lines imprinted on the 
spectra of GRBs or quasars are the main sources of informa¬ 
tion about the chemical and physical properties of high-redshift 
universe. But the bright QSO number is very low at z > 6 
( Fan et al. 20061) . Meanwhile, there are several high-redshift 
GRBs: GRB 050904 at z = 6.29, GRB 080913 at z = 6.7, 
GRB 090423 at z = 8.3 and GRB 090429B at z = 9.4. The 
progenitors of long GRBs are thought to be massive stars, 
so the number of high-redshift GRBs does not decrease sig¬ 
nificantly. The density, temperature, kinematics and chem 
ical abundances can be extracte d from absorpti on lines (Oh 


2002[ Furlanetto and Loeh 2003). For instance, Kawai et al 
(120061) have identified several metal absorption lines in the af¬ 
terglow spectrum of GRB 050904 and found that this GRB 
occur in metal-enriched regions. Two absorption lines have 
been observed in the spectrum of GRB 090423 at z = 8.2 


Figure 19: Observed flux at v = 1.36x 10 l4 Hz (K band) as a function of redshift 
at different observed times, as labelled. The K-band sensitivity of the NIRSpec 
instrument on board the JWST is shown as a horizontal line. The sharp cut¬ 
off at z - 17 is due to Lya absorption in the IGM. (Adapted from Figure 6 in 
I Wang et~ai . (2012ft .) 


(Salvaterra et al. 2009). These lines are due mainly to absorp¬ 
tion metal elements in low ionization stages. 


Wang et al. (2012) studied the ability of metal absorption 


lines in the spectra of Pop III GRBs to probe the pre-galactic 
metal enrichment. The first galaxy simulation carried out 
by Greif et al. ( 2010h was used. The simulation allowed one 
Pop III progenitor star to explode as an energetic supernova, 
then the IGM was polluted by the ejected metals. The simula¬ 
tion box size is 1 Mpc (comoving), and is initialized at z = 99 
according to the ACDM model with parameters: Q m = 1 -Q\ 
Q b = 0.04, h = Hq/ (lOO km s -1 Mpc -1 ) = 0.7, spectral in¬ 
dex n s = 1.0, and normalization erg = 0.9 ( Spergel et al.ll2003l) . 

In Figure [20] the hydrogen number density and metallicity 
averaged along the line of sight are shown within the central ^ 
100 kpc closer to the virialization of the first galaxy at z= 16.4. 
The distribution of metals produced by the first SN explosion is 
highly inhomogeneous, and the metallicity can reach up to Z ~ 
10 -2 5 Z o , which is already larger than the critical metallicity, 
Zcrit < 10 4 Z 0 . Therefore, both Pop III and Pop I/I I stars will 


form during the assem bly of the first galaxies ( Johnson et alJ 


2008tlMaio et~aPl2010h . so simultaneous occurrence of Pop III 


and normal GRBs at a given redshift (Bromm and Loeb 2006; 
de Souza et aklEoTl ). We consider a Pop III burst exploding in 


one of the (still metal-free) first galaxy progenitor minihalos at 
z - 16.4. 

For simplicity, we consider that prior to the GRB only one 
nearby SN exploded beforehand, dispersing its heavy elements 
into the pristine IGM. T wo nucleosynthetic metal y ields for 
Type II core-collapse SNe ( Wooslev and/\Wayerj |F)95 ). a nd for 
pair-i nstability supernovae (PISNe; iHeger and Woos lev 2002, 
20101 ) are considered. Because the hydrogen is substantially 


neutral, metals will reside in states typical of C II, O I, Si II, 
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Figure 20: Possible explosion sites for high-redshift GRBs. Shown are the hydrogen number density (left panel) and metallicity contours (right pan el) averaged 
along the line of sight at z ~ 16.37, when the first galaxy forms. The topology of metal enrichment is highly inhomogeneous. (Adapted from Figure 3 in lWang et al] 


120121 ) .) 


and Fe II, because high-energy photons able to further ionize 
these elements will be absorbed by H I ( Furlanetto and Loebl 
I2003I) . 

Figure|2T|shows two spectra of afterglow at the reverse shock 
crossing time. Top panel is for the top-heavy (Very Massive 
Star) initial mass function (PISN case) and bottom for normal 
initial mass function (Type II SNe case). The cutoff is due to 
Lyman-o' absorption in the IGM which is expected to be still 
completely neutral at z > 10. In the two cases, the metal lines 
are markedly different. The metal yields could be obtained from 
metal lines. So the initial mass function of Pop III stars can be 
derived from the metal absorption lines. 

Strong absorption lines detected in GRB spectra, called 
damped Lyman-a (DLA ) sys tems, could be used to probe 
the metal rnrichment (Savaglio 20061) . For example, 


Castro-Tirado et al. (2013) detected the DLA systems of GRB 


130606A at z = 5.91. The metallicity of this GRB envi¬ 
ronment is in the range from ~ 1/7 to 1/60 of solar from 
the analysis of metal absorption lines. Figure [22] shows the 
metallicities derived from GRB-DLAs and QSO-DLAs. The 
GRB 130606A sub-DLA is only the third GRB absorber with 
sub-DLA HI column density. So GRB sub-DLA is rare and 
hard to find. Simcoe et al. ( 20121) discovered the DLA of ULAS 
J1120+0641 at z ~ 7 with low metallicity. So GRB events at 
Z > 10 from future mission, offer an exciting new window to 
probe pre-galactic metal enrichment in high-redshift galaxies. 


4.3. Cosmic reionization 

The Gunn-Peterson (GP) test tells us that the IGM is almost 
fully ionized at z < 5 ( Gunn and Peterson! 19651) . Reionization 
of the IGM is thought to have occurred during z ~ 6 - 20 by 



Figure 21: Pop III GRB spectrum observed at the reverse shock crossing time 
2© = 16.7 x (1 + 16.4) s. Metal absorption lines are imprinted according to the 
Pop III SN event, PISN vs. core-collapse. The former originates from a very 
massive star (VMS) progenitor, whereas the latter from a less massive one. In 
each case, the cutoff at short wavelengt hs is due t o Lyman -# scattering in the 
neutral IGM. Adapted from Figure 10 in lWang et alj f2012l) .) 


16 


























































o 


1.0 


ffi 

\ 

M CO 




0 2 4 6 

Absorption redshift, z 

Figure 22: The metallicity ([X/H]) as a fu nction of redshift i s shown for 
QSO-DLAs (circles) and GRB-DLAs (stars ISchadv et aj]l201ll : [Thdne et alJ 
120131). including GR B 130606A at z = 5.91 and ULAS J1120+0641 at z ~ 
7 ISimcoe et al.EoT2h . Blue colors are derived for log N( HI) < 20.3 and red is 
derive d for log N(H I) > 20.3. (Adapted from Figure 8 in Castro-Tirad o et alJ 
120131) . 
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Pop III stars and/or quasars, and the precise measurement of the 
reionization i s one of the key topics in modern cosmology (fo r 


reviews, see 


top: 


Barkana and Loe b 2 0011: Robertson et al. 


, 20101 ) . 


The absorption of GRB afterglow is dependent on the struc¬ 
ture of reionization and the global history of reionization, so 
it has the potential to di stinguish between different theoretical 
mode ls of reionization (Miralda- Escudel 1998 : McOuinn et al 


2008). The afterglows at wavelengths close to the Lya res¬ 


onance potentially p rovide a sensitive probe of the ionization 
fraction in the IGM (tMiralda-Escudel7998 ; Barkana and Loebl 
20041 ). The IGM neutral fraction could be derived by fitting 


red damping wing with high precision. The absorption from 
host galaxy complicates the measurement of the IGM ioniza¬ 
tion state from a GRB spectrum, but in principle this absorption 
is less extended in wavelength and could be separated. Some 
studies show that 20-30% of GRB host galaxies have small HI 
column density to allow determination of the absorption from a 
partially ionized IGM ( Chen et ail2007 ). 

By studying the afterglow spectrum of GRB 050904 at z = 
6.3, Totani et al. ( 20061) found that the IGM was already largely 
ionized at z = 6.3, and the upper limit of 0.17 for the neutral 
fraction of IGM at 68% confidence level. But the absorption 
from host galaxy dominates the absorption redward of the Lya 


forest, which limits the constraints on reionization (Totan i et al 


2006; McOuin n et alJl20 08). The bright optical afterglow of 


GRB 130606A at z = 5.9 gives an opportunity to probe the 
ionization status of IGM. The neutral fraction of IGM is found 
to be 0.1 to 0.5 by analyzing of the red Lyo_damping wing of 
the afterglow spectrum taken by Subaru ( Totani et all2013i) . 

From theoretical view, the reionization process can also be 
studied through theoretical model. The average evolution of 
(/ini = n e /nH is derived by numerical integration of the rate 
of ionizing photons minus the rate of radiative recombinations 
(Mad au et al. 1999t Barkana and Loeb 2001 ; Wvithe and Loeb 


N =4000 f >tc =0.1 
N =4000 f asc =0.2 
N =4000 f =0.3 



Figure 23: The HII filling factor <2hii as a function of r edshif t computed for 
different values of / esc - (Adapted from Figure 5 in lWand i2013l) .) 


20031 1 Yu et ap2012h 


^ - a B Cn H Q H n) ^ , (39) 

dz \ n H J dz 

where 

Mon = (1 + z) 3 p*(z)M y /esc/m p (40) 

is the rate of ionizing photons ejected into the IGM, N y is the 
number of ionizing photons, p,(z) is the SFR and / esc is the 
escape fraction. Using the SFR derived from GRBs in section 
13.11 the evolution of the HII volume filling factor (/ H n can be 
numerically calculated from equation ( l39l >. Figure[23]shows the 
evolution of (2 hii as a function of redshift. 

The cosmic microwave background (CMB) optical depth 
back to redshift z is also seriously depend on the reionization 
history, which can be written as the integral of n e cr T d{, i.e.. 


T e (z) = n e (z)o- T ( 1 
Jo 


zT' [c/H(zr>] dz' . (41) 


The optical depth is shown in Figure[24] The WMAP nine-year 
data gives r e = 0.089 ± 0.014 ( Hinshaw et al. 2013 ). which is 
shown as the shaded region. The combination of Planck and 


WMAP data also gives r e = 0.089/® (Planck Collaboration 


2014). So GRB-inferred SFR can reproduce the CMB optical 
depth. But the value of the escape fraction / esc (iRobertson et al 


2Q10i) and clumping factor C are hard to determined. 


5. Summary and future prospect 

GRBs are observed throughout the whole electromagnetic 
spectrum, from radio waves to y-rays, which have been ob¬ 
served in distant universe. Recently, GRBs have attracted a lot 
of attention as promising standardizable candles to construct the 
Hubble diagram to high redshift, as complementarity to other 
cosmological probes, such as SNe la, CMB and BAO. How¬ 
ever, a lot of work is needed to be sure that GRBs can hold 
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Figure 24: The optical depth T e due to the scattering between the ionized gas 
and the CMB photons is shown. The shade region is given by the nine-year 
WMAP measurements. The reionization history calculated from GRB-inferred 
S FR can easily reach r e from WMAP nine-year data. Adapted from Figure 6 
in lWand 1201 3ll .) 


this promise in future. The most important thing is to search 
for a correlation similar to that used to standardize SNe la. 
In order to obtain the correlation, the classification of GRBs 
may be crucial. We must remind that only SNe la are standard 
candles among all SNe. The classical classification method is 
based on the prompt emission properties (duration, hardness, 
and spectral lag). The physics of prompt emission are not fully 
understood ( Zhand 20141), and some new clues from oth er ob¬ 
jects are found ( Wang and Dai 20131 Wang et al. 20141) . But 


sihcation (Zhang 2006; Zhang et al. 2009; 

Lii et al. 2010^_. So 

more physical nature of GRBs is needed 
The circularity problem could be partially 
ing a sample of GRBs within a small reds! 

Zhang et al. 2009). 

solved by analyz- 
lift bin (Lamb et al. 

2005; Ghirlandaet al. 2006a; Liang and Zhang 2006). In par- 

ticular, Liang and Zhang (2006) found tha 

t one can calibrate 

candle correlations 

the power law indices of various standard 


with this method. 


In order to measure high-redshift SFR from GRBs, the re¬ 
lation between long GRB rate and SFR must be known. Be¬ 
sides, theoretical models of the SFR have several free param¬ 
eters, such as the efficiency of star formation and the chemi¬ 
cal feedback strength. From the theoretical SFR, the predicted 
GRB redshift distribution can be derived. So one can use the 
GRB redshift distribution observed by Swift (or future missions 
such as SVOM and EXIST), to calibrate the free parameters. 
More GRB red damping wing with low HI column density are 
required to study properties of IGM. 

Metal absorption lines in the GRB afterglow spectrum, giv¬ 
ing rise to EWs of a few tens of A, which may allow us to 
distinguish whether the first heavy elements were produced in a 
Pop III star died as a PISN or a core-collapse SN. To this extent, 
the spectrum needs to be obtained sufficiently early, within the 


first few hours after the trigger. Upcoming JWST would detect 
much more high-redshift GRBs (properly Pop III GRBs) with 
high resolution NIR spectra including metal absorption lines, 
which allow one to measure the cosmic metallicity evolution. 

In the future, the French-Chinese satellite Space-based multi¬ 
band astronomical Variable Objects Monitor (SVOM) and 
JWST, have been optimized to increase the number of GRB 
and the synergy with the ground-based facilities. There are a 
combination of multi-wavelength detectors on board of SVOM 
(Paul et al. 2011). ECLAIRs wide-held camera will detect 
GRBs in the energy range of 4-150 keV. The spectral infor¬ 
mation of prompt emission will be measured by Gamma-Ray 
Monitor (GRM). The afterglow can be obtained by the Micro 
channel X-ray Telescope (MXT; 0.3-10 keV) and the Visible 
Telescope (VT; 400-900nm). SVOM can detect about 80 GRBs 
per year, and more than 50% of GRBs have redshift measure¬ 
ment (Petitjean and Vergani 2011). JWST is a large, infrared- 
optimized space telescope with 6.6 m diameter aperture. It 
has four scientific instruments: a Near-IR Camera (NIRCam), 
a Near-IR Spectrograph (NIRSpec), a near-IR Tunable Filter 


Imager (TFI), and a Mid-IR Instrument (MIRI) (Gardner et al 


2006). But the direct detection of a single Pop III star is not fea¬ 


sible even for JWST, i.e., the AB magnitude of a M - 1000M o 
star is only 36 at z ~ 3 0. Meanwhile, the Pop III GRBs can 
be detectable by JWST ( Wang et all 2012k Mesler et al. 2014 ; 


Mac pherson et al.ll20 13). This will boost the amount of infor¬ 


mation available to tackle the important issues revealed by this 
exciting held of research. 
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